Correlation, excitation, and relaxation of the antiferromagnetic nanoparticle 
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We study the dynamics of spins in the canted antiferromagnetic hematite nanoparticles a-Fe2 03. 
The system includes an antiferromagnetic exchange (J) between two sublattices as well as the 
potential dominated by a uniaxial anisotropy and dissipates through the thermally triggered spin- 
phonon mechanism. The exchange J is found to reduce the superparamagnetic relaxation and drive 
the incoherent relaxation of transverse spins in the low temperature. Dynamical properties are semi- 
quantitatively understood within a role of the exchange interaction J and agree well with recent 
inelastic neutron scattering experiments. 



Magnetic properties of the fine particles of a size ~ 
O(10) nm are interesting and currently attracting much 
attention. In the system, the most relevant energy scale 
tends to be the magnetic anisotropic energy, easily com- 
parable to the thermal energy T, so that the magneti- 
zation (or superspin) is subject to considerable thermal 
fluctuation, called the superparamagnetic relaxation^ . 
Phenomena of such magnetic fluctuations have the po- 
tential importance in the fundamental research as well 
as the technological application 2] . 

Hematite, a-Fe203, is a mineral of the corundum 
(rhombohedral) structure at room temperature, charac- 
terized by the rhombohedral axis [111] and the basal 
(111) plane. Spins in hematite are sightly canted away 
from antiferromagnetic (AFM) orientation on the basal 
plane, giving place to the weakly ferromagnetic (FM) 
state, between the Morin temperature Tm ~ 260 K and 
the Neel temperature Tjy ~ 956 K, but spins are ori- 
ented along the [111] axis and perfectly antiferromag- 
netic below Tm- Magnetic properties are found to depend 
on the particle size and Tm is less than 5 K in typical 
nanocrystalline hematites 3] , where several kinds of mag- 
netic anisotropy energies, not present or much smaller 
in the bulk, play important roles. In nanocrystalline 
hematites, the magnetization oscillates due to superpara- 
magnetic relaxation along the easy axis and collective 
magnetic excitations in vicinity of the easy axis[2j. 

Many studies of the magnetic fluctuations in the 
nanoparticle hematite are carried out using Mossbauer 
spectroscopy or measuring the magnetic susceptibility. 
Morup and his coworkers have used Mossbauer spec- 
troscopy and studied the energy scale of the superpara- 
magnetic relaxation and anisotropic barrier 0. They 
have also explored effects of the inter-particle interac- 
tion in strongly interacting nanoparticles (e.g., uncoated 
dry particles or a form of mixture such as a-Fe 2 03+CoO 
or 57 Fe-NiO+CoO) and found the absence of fast super- 
paramagnetic relaxation, i.e. the slowing down of the 
relaxation^. By the way, the magnetic relaxations in 
the system are argued to occur highly fast, being almost 
the ultrafast process (in the time scale of < 0(1O~ 12 ) 
seconds or in the frequency scale of THz range), and the 
Mossbauer spectroscopy or susceptibility measurements 
are then found to be slow probes for the relevant mag- 
netic dynamics. By the same group and coworkers, due 



to those reasons, the inelastic neutron scattering experi- 
ment is suggested as one of valuable complements @- 

In this paper, we investigate the detailed dynamics of 
the magnetization in the AFM nanoparticle system. The 
system incorporates the AFM exchange J between two 
sublattices as well as the potential by the anisotropy and 
undergoes the spin-relaxation through the spin-phonon 
mechanism. A main goal of the paper is to understand 
the respective dynamics and relaxation of both the su- 
perparamagnetic fluctuations and the collective magnetic 
excitations with respect to the temperature T. Prior to 
the present work, theoretical studies have analytically 
shown that the superparamagnetic relaxation follows the 
simple Arrhenius behavior of ~ To(T)e~ EB / T with a cer- 
tain form of Tq(T), where Eb the energy barrier, for 
FM grains (like Mni 2 0i 2 molecules)0 H 0. At the 
very low temperatures (T < 2 K), however, it is found 
that the relaxation deviates from the Arrhenius law and 
is dominated by a resonant tunneling between low-lying 
states 10]. In the present work, such low temperature 
regimes are not assumed. 

Considering that spins in the bulk hematite are in the 
slightly canted AFM (or weakly FM) state on the (111) 
basal plane, where the hexagonal crystal anisotropy is 
dominated by the uniaxial one, we suggest the Hamilto- 
nian TIq as 



Ti-o = K Z (S 2 Z 



Si) 



Ky(S 2 



S 2y) 



JSr-Sa, (1) 



where Si and S 2 are the magnetization or superspin 
(or simply spin) of two AFM sublattices, with a size of 
S>1, and K z and K y are uniaxial anisotropic constants 
along z ([111] direction) and in the basal plane, respec- 
tively, and J is the exchange between two sublattice mag- 
netizations. K z is expected to be the largest energy scale, 
i.e. K z ^> K y , and K z ^> J. K Z S 2 is noted to be a natu- 
ral energy unit. Hq can be easily diagonalized by taking 
the basis of \S\ Z ,S2 Z ) and performing the calculation of 
Clebsch-Gordan coefficients for the coupling term. In ad- 
dition to H.q, we consider the spin-phonon interaction V 
(= Vi + V 2 ) raising the spin- relaxation, Vi — Vif, 



v i — 9o{Si z Si x + SiySi, 



) + giS ix Siy + h.c. 



(2) 



where / denotes the lattice deformation. The first term 
in Eq.@ gives rise to transitions from \rrii) to \m,i ± 1), 
while the second term leads to transitions to \mi ± 2). 
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FIG. 1: Left panel: respective spin averages divided by S(S + 
1) with respect to T. £ = 0.3 and ( — 0.5. Right panel: 
behaviors of ((Si - S 2 ) 2 } - S? - S| normalized by 2S(25 + 1) 
with respect to T. 



We have interests in the relaxational dynamics of spins 
at the AFM Bragg scattering, which is given by the fol- 
lowing correlation function, = Si^ — S^, /x = x, y, z, 



-iCr 



Sn 



(3) 



where C is the Liouville operator, CO = [H,0], 
H = Ho + V. It is the spectral function J PM (w) — 
J dre iwr (/) MAl (r) that can be compared directly with the 
inelastic neutron scattering experiment. Up to the sec- 
ond order of V, the correlation function (/^(t) is written 
as ^(t) = (e m ° r S' A1 e- mor e- 4£vT S' Al ), which is valid to 
the lowest order in the memory function|ll|. By intro- 
ducing the relaxation matrix 7Z whose elements consist 
of the transition probabilities of the Markov equation, wc 
replace e~ lCvT by e~ UT 



Mr) = (e m ° T S^ 



-Ut 



(4) 



With the basis of the eigenstate of 7io, elements of 
the relaxation matrix TZ are found as follows; 7Z nn — 
Eln+i 7n-m + ELo e _(£ "'^" )/T 7m-.n, and 1Z nm = 
-7m->n for n > to, and 7Z nm = - e -(- E «- B ™)/' 7 ^ n _ m for 
n < to, where W(= (25 + l) 2 ) is the dimension of the 
matrix and finally 7„_> TO is defined by@, 



7n 



|(n|i;|m)|'' 



1 _ e -(E n -E m )/T 



(E n — E m ) 3 



We assume E n > E m for n < to. TZ is in fact 
a generalization of Orbach's process [12j to a spin cas- 
cade. For a coupling to elastic waves, 7„_ >rn should 
be oc (E n — E m ) 3 from the line shape function. It 
should be noted that the eigenstate of the relaxation 
matrix 1Z does not construct an orthonormal basis set 
because 1Z is not a (real) symmetric matrix [T^. We 
then need an orthogonalization process of nonorthogonal 
basis, i.e. like the Gram-Schmidt orthogonalization[l4j. 
Putting the eigenstate of 1Z with the relaxation eigen- 
value of Xi as \pi), the orthogonalized state \tpj) is, \ipj) = 
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FIG. 2: Left panel: total spectra J(w), a sum of J X x{ui), 
Jy y {uj), and J zz (uj) for given parameter sets with respect to 
T. Right panel: transverse spectra J yy (uj) and Jzz(uj) and 
their sum J±(u). The unit of uo is 2K Z S. 



with 



bo), 



i/Vi-EtoKPiki)! 2 . and c t = fOT * > h 

and we find Cj = — Cj E;=i Q (PilVi) f° r * < 3- Now 
the correlation function of (j)^ (r) is expressed as 

Mt) = EEe-^e^-^H^TO) (5) 

n m 

x EE c ] e ~ AlT H^)<^is>>- 

j i=0 

In the left panel of FigQ] it is shown that x should be 
the direction of the easy axis, while, in the right panel, 
the temperature variation of AFM correlation strengths 
are given for £'s and £'s (£ = K y /K z and £ = J/K z ). 
Energy quantities including T are given by a unit of 
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FIG. 3: (a) the excitation energy of S y (A y ). Dashed lines 
are fittings from oc \/£C- (b) the excitation energy of S z (A z ), 
oc \/F+"C The legend in (a) is applied for (b), too. The 
inset shows the transverse excitation energy (Ax = A y = 
A z ) when £ = 1, nicely fitted with oc y/£ + 1. Results are 
evaluated at T — 0.1. The unit of both A y and A z is 2K Z S. 



K Z S 2 , but the excitation energy u> (also A y and A 2 
given later) should be scaled by 2K Z S rather than K Z S 2 
from K Z S 2 - K Z (S - l) 2 w 2if z S*. At low temperatures 
(T <C K Z S 2 ), there are two distinctly different types 
of dynamics, namely the superparamagnetic fluctuation 
J\\(u>) (= J a; x(^-')) and the collective magnetic excitation 
J±(oj) (= J ra (u>) + J zz (uj)). 

Collective magnetic excitations lead to the inelastic 
scattering process. The relative intensity of J yy (uj) 
and J zz (lu) depends on the ratio of two different spin- 
dissipation channels in Eq.JSJ. We take the ratio gi/go = 
0.5 because transitions of Srrii = 1 are considered to oc- 
cur more probably than Srrii — 2 and ag 2 = 1 x 10~ 3 . 
Under the ratio, it is found that J±(u>) is governed by 
J yy (uj), J±(u>) ~ J yy (uj). For three sets of parameters, 
we provide the spectral functions in FigEl We note that 
two important features characteristically develop with re- 
spect to £, £, and T; one is the collective excitation en- 
ergy and the other is the relaxation of spins. By com- 
paring panels in Fig|2 it is found that the excitation 
energies of S y (A y ) and S z (A z ) behave differently. In 
FigI3 the collective excitation energies for £ and £ are 



Aj_) is scaled roughly by 
/£C for £ < C and A y oc £ 



given. We find that A y (? 
~ that is, A y cx 

for ( ~ 0. This is consistent with the result of inelas- 
tic mode energy of a pure antiferromagnetj^, Il5| . On 
the other hand, A z is by ~ VC+T- Especially, by the 
symmetry, for a case of £ = 1, A y should be same as 
A z , Aj_ - A y = A z cx VC^' 



A y goes to by J yy {uj) = J. 



F 1, 



while, for £ 



0, 



Now, let us con- 



sider the relaxation of the longitudinal and transverse 
components of spins, which is in principle understood 
from the broadness of curves in Fig0 Instead of di- 
rectly measuring the broadness, however, for the system- 
atic investigation of the relaxation, we define and pro- 
pose a useful renormalized dissipation parameter 
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FIG. 4: Behaviors of spin relaxation constant fx(T), y y (T), 
and 7z(T) for £'s and f's with respect to T. Figures (a) and 
(b) are for £ = 0.5 and (c) and (d) are for £ = 0.5. Thick 
solid lines in (a) and (b) are for £ = 0. The legend in (a) is 
also for (b) and that in (c) also for (d), respectively, (e) gives 
the energy barrier out of the Arrhenius law y x (T) ~ e ~ E Bl T 
with respect to £ (the line is a guide for the eye). 



where 7^ m is found to describe the relaxation of the 
system at least in an approximate sense. It is evident 
that the relaxation of the system is mainly governed by 
some smallest eigenvalues of A^'s and it is then essential 
to account for the long-time behavior after long enough 
r has elapsed. Due to the reason, we obtain j£ m by 
integrating both sides from r = to t = oo, j£ m « 

(m|S»/ [E j Eio {C i i WPi)(<Pi\SM/fri + S )}] ~ S > 
where the finite dissipation 8 = 2.5 x 10~ 4 K Z S 2 has 
been added for the stable calculation. 8 can be 
understood corresponding to the finite energy reso- 
lution in an actual experiment 16] . With 7# m , the 
time-dependent spectral function MM (r) is MM (r) = 



y t 

Z—m Z— /r. 



-0E ne i(E n 



- E ™+ l < m ) T \(n\S Jm>| 2 , from which 



7 M is further defined as the smallest value of 7^ m at 
the peak position of the spectral function of 5 M . Be- 
haviors of 7 M (T) have essential importances in that the 
magnetic relaxation of spins and related dynamics of en- 
ergy dissipation can be well understood through them. 
We have two distinguishable kinds of the spin relaxation, 
as shown in Fig^J one is the superparamagnetic relax- 
ation along the easy axis (j x (T); FigsQJa) and (c)) and 
the other is the relaxation by the collective excitations 
transverse to the easy axis (■y y (T) or j z (T); Figs^Jb) 
and (d)). Because j z (T) is always larger than 7 y (T), 
consistent with J yy (u>) 3> J zz (uj), it is reasonable to put 
7_i_(T) = J y (T). Like cases of FM grains, the superpara- 
magnetic relaxation still follows the simple Arrhenius be- 
havior as a first approximation, even if the relaxation has 
two-dimensional nature. However, the energy barrier de- 
pends on £ and £ and does not have so simple a form. 
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From Fig^Ja), it is evident that the AFM exchange J re- 
duces the superparamagnetic relaxation, in other words, 
the energy barrier increases with J as shown in FigBJe). 
This finding is consistent with the claim of the absence of 
the fast superparamagnetic relaxation in the interacting 
nanoparticles through the Mosbauer spectroscopy. With- 
out the exchange, the energy barrier approximately fol- 
lows ~ K y S 2 , which implies ^ x cx e~^ T ; j x cx e ~°- 46 / T 
for £ = 0.5 and £ = 0. Figures Efb) an d (d) give the 
transverse relaxation, quite different from the superpara- 
magnetic relaxation. In the low temperatures, to increase 
the temperature reduces 7_i_(T). 7_i_(T) has then a min- 
imum and slowly increases with T. In such a sense, the 
transverse responses of spins are found incoherent in the 
low temperatures. We conclude that the incoherence is 
driven by the AFM spin exchange J because 7_l(T) is 
always coherent for J = (a thick line in Fig^b)). The 
incoherent behaviors of ~f±(T) in the low temperatures 
has been actually observed in the inelastic neutron scat- 
tering experiment for a-Fe203@|. 




Their estimations have been done with appreciable er- 
ror ranges from temperature fits in various ways by 
Morup and his coworkers 6] . A typical one of them is 
K y S 2 ~ 500 K (i.e. ~ 42 meV) from Arrhenius behavior 
of the superparamagnetic relaxation, but it is found from 
FigEIc) that the value may have been overestimated al- 



most by 60 %. Thus, we estimate K y S 2 ~ 300 K (~ 25 
meV), while the value of K Z S 2 is more uncertain because 
its value is not known for the nanoscopic system, but of a 
same order of magnitude with the bulk case. So we take 
K Z S 2 ~ 1000 K (~ 83 meV), that is, we have £ = 0.3. On 
the other hand, the unit excitation energy 2K Z S can be 
found ^0.5 meV to be consistent with the experiment. 
The inelastic energy mode (Aj_) is then ~ 0.3 meV at 
T = 0.1 (~ 100 K). The value of JS 2 can be estimated 
as 300-500 K from the transverse excitation energy Aj_ 
in Fig|3| and from the relaxational behaviors of spins in 
Fig0] Based on the estimations, J X (T) is about O(10 -3 ) 
meV and 7_l(T) about O(10 _1 ) meV in the temperatures 
T < 0.2 (i.e. T < 200 K). We note jx(T) is larger than 
7a:(T) by O(10 2 ), which is consistent with the actual sit- 
uation of the experiment 6] . 

In summary, we have investigated the dynamics of 
spins in the antiferromagnetic hematite nanoparticle a- 
Fe2C>3 with the explicit AFM exchange J. The exchange 
J has been found to reduce the superparamagnetic re- 
laxation and also induce the incoherent oscillation of 
transverse spins in the low temperature. The spin re- 
sponses are well understood within the model and a semi- 
quantitative estimation of the relevant parameters has 
been done through a comparison with the recent inelas- 
tic neutron scattering experiments. 
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